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ABSTRACT 

We report the application of the new Monte Carlo method, Smoothed Particle Inference (SPI, de- 
scribed in a pair of companion papers), towards analysis and interpretation of X-ray observations 
of clusters of galaxies with the XMM- Newton satellite. Our sample consists of publicly available 
well-exposed observations of clusters at redshifts z > 0.069, totaling 101 objects. We determine the 
luminosity and temperature structure of the X-ray emitting gas, with the goal to quantify the scatter 
and the evolution of the Lx - T relation, as well as to investigate the dependence on cluster substruc- 
ture with redshift. This work is important for the establishment of the potential robustness of mass 
estimates from X-ray data which in turn is essential towards the use of clusters for measurements 
of cosmological parameters. We use the luminosity and temperature maps derived via the SPI tech- 
nique to determine the presence of cooling cores, via measurements of luminosity and temperature 
contrast. The Lx — T relation is investigated, and we confirm that Lx oc T^. We find a weak redshift 
dependence (ex (1 + z)^^'^ ^ Plt = 0.50 ± 0.34), in contrast to some Chandra results. The level of 
dynamical activity is established using the "power ratios" method, and we compare our results to 
previous application of this method to Chandra data for clusters. We find signs of evolution in the 
P3/P0 power ratio. A new method, the "temperature two-point correlation function," is proposed. 
This method is used to determine the "power spectrum" of temperature fluctuations in the X-ray 
emitting gas as a function of spatial scale. We show how this method can be fruitfully used to identify 
cooling core clusters as well as those with disturbed structures, presumably due to on-going or recent 
merger activity. 

Subject headings: galaxies: clusters: general — X-rays: galaxies: clusters 



1. INTRODUCTION 

Clusters of galaxies are the largest gravitationally 
bound structures in the Universe, and thus should pro- 
vide a fair sample of its matter content. This makes clus- 
ters good candidates for cosmological studies. In partic- 
ular, the gravitational growth of initial density pertur- 
bations can be used to constrain cosmological parame- 
ters via determination of the mass function of clusters of 
galaxies (e.g. Voit 2005), but this requires good knowl- 
edge of cluster masses. 

One of the most promising avenues towards the mea- 
surement of the mass function of clusters is based on es- 
timates of the cluster temperature and luminosity from 
X-ray observations of a large number of objects. Specifi- 
cally, the number density for clusters of different masses 
can be estimated using mass-observable relations, cal- 
ibrated using nearby clusters where spatially resolved 
spectros copy is a vailable. The mass-temperature (M- 
T) (e.g. lArnaud et al. 20 05|) and lumin osity-temperature 



{Lx-T) (e.g. lArnaud fc EvrardI 11999 ^ scaling relations 



are of particular importance since temperature is often 
used as a proxy for mass and the relation to luminosity is 
needed to understand the sample selection function since 
X-ray selected samples generally are flux limited. More 
recently, a new proxy for cluster mass, Yx, has been pro- 
Electronic address: kanderss@space.mit.edu 



posed (|Kravtsov et al.l [20061 ). This quantity is a simple 
product of the gas mass of the intra cluster medium and 
its temperature and has shown to exhibit a low amount 
of scatter. 

Clusters are formed hierarchically through mergers of 
smaller clusters and groups. This merging activity is 
observed as distortions of their X-ray surface bright- 
ness profiles. Substructures and mergers affect the mass 
determinations and increase the scatter in the scaling 
relations. Measurements us ing only "relaxed" clu sters 
have achieved high precision (|Vikhlinin et al.|[2006af ) and 
selective studies, such as those of the gas mass frac- 
tion (fgas) in clusters focus only on t he largest cluster s 
with minimal amount of substructure ([Alien et al.ll2008l ) . 
However, substructure was found to be present in ^ 50% 
of clusters in a ROSAT study (jSchuecker et all 120011 ) 
and studies of XMM data in a REFLEX-DXL study 
find substructures prese nt in all clusters in that sample 
(|Finoguenov et ani2005l ). Understanding, and assessing 
the effect of substructure on the robustness of mass de- 
termination is thus crucial. 

Besides the complications associated with the merger 
activity, evolution of clusters can affect the applicabil- 
ity of scaling relations. The high density environments 
in galaxy cluster cores cause them to cool radiatively 
and this is observed in undisturbed clusters as a decrease 
of the average projected temperature towards the center 
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(e.g. iPeterson fc Fabianl[2QQ6h . The undisturbed "cool- 
ing clusters" have sharply peaked luminosity profiles, a 
feature that is not observed in their non-cooling counter- 
parts. Cool-core clusters deviate from the Lx-T relation 
since the core has higher luminosity and lower temper- 
ature than the cluster population on average. Obser- 
vational evidence exists supporting the argument that 
cool cores can survive to some extent during a clus- 
ter merger. These core remnants are then observed 
as sharp contact disc ontinuities or "cold fronts" (e.g. 
iMarkevitch fc Vikhlinin 2007), possibly affecting the M- 
T scaling. 

All this indicates that it is important to assess and 
quantify the dynamical state of a cluster when the cluster 
data are used for the determination of cosmological pa- 
rameters. While this has been success fully attemp ted in 
the past (see e.g lMaughan et al.|[2QQ8l : iChen et al.ll2007[ ). 
statistical studies of cluster substructure have so far only 
used the information from the spatial distribution of X- 
ray counts on the sky, hence only mapping the luminosity 
structure. Ideally, the knowledge about the temperature 
structure should also be included in these searches since 
it holds important information about the dynamical his- 
tory of the cluster. 

The Smoothed Particle Inference (SPI) tech- 

nique, developed recently by some of us 

([Peterson. Marshall, fc AnderssonI l2QQ7l ) as an al- 
ternative to standard analysis techniques, is well suited 
for detecting the effects of cooling cores and substructure 
from both the projected temperature and luminosity 
distributions of clusters (see Andersson et^aL 2007) . 
This method relies on a description of a cluster as a 
large set of smoothed particles (two-dimensional, spatial 
Gaussians), each of which is described by a luminosity, 
spatial position, Gaussian width, temperature, redshift, 
and a set of elemental abundances. A large set of these 
particles is propagated through an instrument model, 
and the model parameters are adjusted using Markov 
Chain methods. The resulting distribution becomes a 
kilo-parametric description of the cluster. 

In this paper, the SPI method is used on a large num- 
ber of cluster observations available through the XMM 
public archives. A cluster model is built using the the 
imaging spectroscopy XMM data. The output of the 
modeling is used to separate clusters with cooling cores 
from more disturbed clusters and to study their prop- 
erties separately. Specifically, we aim to assess the ef- 
fects of cooling cores and substructure on the luminosity- 
temperature relation and study any possible redshift de- 
pendence of these effects. We also apply a new statistic, 
including the spatial distribution of both luminosity and 
temperature to quantify the level of dynamical instability 
present in the clusters. This statistic is designed to dis- 
tinguish cool core clusters and isothermal clusters from 
those with more temperature structure indicating a re- 
cent or ongoing merger event. 

In Section 2 we describe the construction and proper- 
ties of the cluster sample and outline the data reduction 
scheme, in Section 3 we explain the different methods 
used to analyze the data, and in Section 4 we describe 
the processing performed on the output cluster models 
and the methods used to quantify cluster properties. In 
Section 5 we display the statistical results of our model- 
ing and in Section 6 we discuss the possible systematic 



effects associated with this. We finally conclude in Sec- 
tion 7 with a discussion on the results, problems and 
possible improvements. 

In all calculations we have assumed a concordance 
cosmology with Qm = 0.3, = 0.7 and Hq = 

70 km s-^ Mpc-^ 

2. CLUSTER SAMPLE AND DATA PREPARATION 

The cluster sample was compiled by cross-correlating 
the NASA Extragalactic Database (NED) with the pub- 
lic XMM-Newton observation archive as of May 10, 2006. 
The requirement for selection in NED was that the clus- 
ter should be a known X-ray source and a known galaxy 
cluster (2005 clusters) or group (120 groups). The XMM 
pointing was required to be within 3.5' of the source po- 
sition in NED. This resulted in 278 matches for clusters 
and 34 for groups, some of which were multiple matches 
of the same cluster. After multiples and sources not vis- 
ible in the X-ray data were removed the total number of 
clusters and groups was 201. We further removed clusters 
where we could not get a reliable spectrum, i.e. clusters 
with a fluence (time-integrated energy flux) below 10~^ 
erg cm~^, as well as nearby sources where we could not 
fit 1 Mpc within a 13' radius. This selection left us with 
101 sources (see Table [T]). We note that although this 
sample is by no means complete, it represents a broad 
range of cluster properties over a large range of redshift. 

The sample encompasses clusters with luminosities 
from 9 X 10^^ erg s~^ to 1.1 x 10^^ erg s~^, average tem- 
peratures from 2.2 to 11.6 keV (See Section 13. ip and 
redshifts from z = 0.069 to z = 0.89. The distribution 
of luminosities with redshift for the sample is shown in 
Figure [T] where the names of a subsample of well-known 
clusters are printed in the plot. Luminosities were cal- 
culated from the observed flux, using tabulated redshifts 
(as listed in NED), assuming a concordance cosmology 
(mentioned above), and applying a bolometric correc- 
tion. 

2.1. Data reduction 

The data were reduced using standard pipeline process- 
ing and the calibration implementation as of XMM Sci- 
ence Analysis Software (SAS) version 6.5. Background 
flares from soft protons were removed using lightcurve 
filtering in both soft (MOS: 0.3-10 keV, PN: 0.3-12 keV) 
and hard (MOS: 10-12 keV, PN: 12-14 keV) X-ray bands. 
For the soft band, light curves were binned in 10 s inter- 
vals while we used 100 s bins for the hard band. In both 
cases the data were not included during the time when 
the count rate exceeded 3<j above the quiescent count 
rate, indicating a proton flare. 

The event files were also filtered for non- X-ray events 
by selecting only single and double pixel events for PN 
and single to quadruple events for MOS. Bad pixels 
and pixel columns were removed by applying the stan- 
dard keywords in event selection; FLAG=0 and #XM- 
MEA_EM / #X MMEA_EP. The data re duction follows 
that described in I Andersson et al] (|2007[ ). 

3. DATA ANALYSIS METHODS 

We analyze all objects using both simple spectral anal- 
ysis, where a single spectrum is extracted and fitted, as 
well as the SPI analysis, using the Monte Carlo approach 
and modeling the clusters both spatially and spectrally. 
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TABLE 1 
Spectral fit parameters 



Name 


z 


T [keV] 


Abundance wrt. Solar 


Lhol [er 


gs 1] 


Exposure [ks] 


CIZAJ1645.4-7334 


0.069 


3.861^22 


0.65l°i? 


3.9 


1044 


14 


A1837 


0.070 


3-85l°-°^ 


0.38l°-°3 


1.7 


1q44 


45 


A3112 


0.070 


/, OO + 0.05 


0.501^2 


6.4 


1q44 


21 


A1775 


0.072 




0.49+°-°3 


2.9 


1044 


23 


A399 


0.072 


6.92+°-23 




7.4 


1044 


11 


A1589 


0.072 


4.8ll°;?g 


0.30l°-°6 


2.1 


1q44 


15 


A2065 


0.073 


5.36_Q -|^Q 




6.4 


1044 


18 


A401 


0.075 


7 79+0:22 

• • ' Q 23 


PI 97+0.03 


1.4 


1045 


11 


A2670 


0.076 






3.3 


1044 


17 


A2029 


0.077 


6.88^Q'Qg 


0-411^2 


2.2 


1045 


11 


RXCJ1236.7-3354 


0.080 


2.8010-lS 


0.46l°-°? 


1.6 


1044 


12 


RXCJ2129. 8-5048 


0.080 


4.55i°;f9 


PI 97+0.06 

yj-^ 1 Q Qg 


2.6 


1044 


21 


A2255 


0.080 


6-98to28 


0.281°-°^ 


6.7 


1044 


10 


RXCJ0821. 8+0112 


0.082 




0.28lg-S° 


1.3 


1044 


9 


RXCJ1302. 8-0230 


0.083 


cro + O.OQ 

O.J^_Q Qg 


0.491^6 


1.9 


1044 


22 


A1650 


0.084 


^ P-o + O.'O? 

— 08 




7.8 


1044 


37 


A1651 


0.084 


6.20+°-" 

— U. ib 


0.34+°03 

— U.Uo 


10.0 1044 


10 


A2597 


0.085 


3.46l°°3 


0.39l°°l 


6.0 


1044 


56 


A1750 


0.086 






3.1 


1044 


29 


A478 


0.088 




0.37l°-°i 


2.7 


1045 


96 


A278 


0.089 


3.39l°;j| 


0.261^^5 


1.5 


1044 


27 


A2142 


0.090 


8.15 Q 2g 




2.9 


1045 


6 


A3921 


0.094 


5.65+°-}! 

— 0.14 


0.33+°°l 

— 0.04 


6.3 


1044 


29 


A13 


0.094 






2.8 


1044 


31 


A3911 


0.097 


5.94+°" 

— U. 14 


0.26+°°? 

— 0.03 


5.7 


1044 


25 


RXCJ2319. 6-7313 


0.097 


2 27+0.09 

— 0.07 


0.32+°°', 

— 0.05 


2.2 


1044 


8 


CL0852+1618 


0.098 




— U. i9 


9.6 


1043 


31 


A3827 


0.098 




27+°-°2 

yJ-^ i Q Q2 


1.1 


1045 


21 


RXCJ0211. 4-4017 


0.101 


9 99 + 0.07 
— 07 


43+0-07 

— 06 


1.2 


1044 


26 


A2241 


0.101 


3.66+°-i 

— U.28 


0.81+°-?« 

— 0.18 


3.0 


1044 


5 


PKS0745-19 


0.103 


6-441^7 


0.34l°°l 


3.4 


1045 


20 


RXCJ0645. 4-5413 


0.105 


1 QQ + 0-33 
' •'-*'^ — 26 


o.24ini 


7.7 


1044 


13 


RXCJ0049. 4-2931 


0.110 


4.02+0-" 

— U. io 


0.40+°°6 


3.4 


1044 


19 


A1302 


0.116 


6.59l°;« 


— 10 

0.53_Q Qg 


5.6 


1044 


16 


RXCJ0616.8-4748 


0.116 


5-05_Q 


0.26l°i° 


3.0 


1044 


9 


RXCJ2149. 1-3041 


0.118 


crq + 0:07 

— 07 


0.47l°-°l 


3.8 


1044 


23 


RXCJ1516.3+0005 


0.118 


9c: + 0:i6 
— U. io 


0.30+Ht 


4.7 


1044 


25 


RXCJ1141.4-1216 


0.119 


p-q+o.oe 


Q-Q^ 

0.53 Q Q2 


4.1 


1044 


26 


RXCJ0020. 7-2542 


0.131 


6.47+°-2^ 

— U.Z6 


0.24+°-°^ 


5.8 


1044 


16 


RXCJ1044.5-0704 


0.134 




0.361°-°^ 


7.3 


1044 


25 


RXCJ0145.0-5300 


0.136 




0.34l°:°« 


8.3 


1044 


14 


A1068 


0.138 




39+°-25 

U.Oi?_Q Q3 


8.4 


1044 


20 


RXJ1416.4+2315 


0.138 


f:o + 0.38 
<J-Oo_0.33 


n 25+0-14 

U.^O_Q 13 


2.2 


1044 


6 


RXCJ0605.8-3518 


0.141 


^•O^_0.08 


QQ + 0-03 

U.O»_Q Q3 


1.0 


1045 


21 


A1413 


0.143 


7 on + 0.19 


36+0-04 

U.oD_Q Q4 


1.4 


1045 


24 


RXCJ2048. 1-1750 


0.147 


5 Q2 + 0.28 
"^•^^-0.23 


p, 90 + 0.05 

U.ZO_Q Q5 


5.8 


1044 


23 


A3888 


0.151 


Q 01 +0.69 

^•"^-^-o.si 


p, 90+0.08 
u.zo_Q Q7 


1.5 


1045 


5 


A2034 


0.151 


7 41+0-27 


29+0-04 

^•^^-0.04 


1.6 


1045 


12 


RXCJ2234.5-3744 


0.151 


7 00 + 0.22 


lQ + 0-03 
'-'•-^^-0.03 


1.5 


1045 


23 


A2204 


0.152 


6 44+0-08 


n 07+0.02 
-0.02 


3.5 


1045 


19 


RXCJ0958.3-1103 


0.153 


9r) + 0.19 


40+0-05 

U.^U_Q Q5 


1.0 


1045 


8 



Note. — Redshifts, temperatures, metal abundances (wrt. Solar), luminosities and exposure times 
for the sample. Redshifts are taken from the NED listing. Temperatures and metal abundances were 
derived using spectral fits (Section I3.1I ). Luminosities within a 1 Mpc aperture were derived using the 
SPI modeling (Section [3T2l . 
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Fig. 1. — The distribution of bolometric luminosities for the sample. Low Lx clusters (< 4 lO^^erg s ^) are shown as red circles, 
intermediate Lx clusters (4 10^^ erg s"-*^ < Lx < 2 10^^ erg s"-*^) as green stars and high Lx clusters (> 2 10"^^ erg s"-*^) as blue squares. 
Simulated cool core clusters are shown as gray filled circles. 
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TABLE 2 
Spectral fit parameters, Contd. 



IN clIIlG 


z 


1 [KG V J 


AbunQcincG wrt. oolcir 


-ll 

^boi [erg s J 


n 1 

Exposure [ksJ 


A868 


0.153 


^•^"-0.39 


"^•^-^-0.08 


6.8 10^4 


6 


RXCJ2014.8-2430 


0.154 


^•^^-0.07 


n Q7+O.02 
-0.02 


2.0 10^^ 


23 


A2104 


0.155 


^•^^-0.97 


-un 1 zL 

0.38lg:l^ 


1.3 10^^ 


5 


RXCJ0547.6-3152 


0.166 




n '^l +U.U4 
"^•'^-'^-0.03 


1.1 10^^ 


22 


A2218 


0.171 


7 1 7'+0.21 
-0.17 


n 99+0.03 


1.1 10^^ 


16 


A1914 


0.171 


, -Un 9R 

^•^^-0.18 


"^•^^-0.03 


2.7 10^^ 


22 


A665 


0.182 


7 04+0-23 
' -^^-0.24 


n 9p;+0-03 
"^•^^-0.03 


1.7 10^^ 


57 


A1689 


0.183 


„ „_-|-0 17 

9.07lg:i^ 


^ ^^-1-0 02 


3.1 10^^ 


35 


A383 


0.187 


4 21 +0-09 
^•^-^-0.08 


4fi+0-04 
^•^"-0.03 


9.7 10^4 


28 


A520 


0.199 


_ . --I-D .3.3 

8.45lg-J^ 


^ _ . _i_n 04 


1.7 10^^ 


38 


A2163 


0.203 


11 10+0.36 

J_J_. J_Z_Q 3g 


91+0-04 

'-'•^-'--0.04 


5.7 10^^ 


10 


A209 


0.206 


7 91 +0-27 
' -^-^-0.26 


n 27+0.04 

-0.04 


1.5 10^^ 


17 


A963 


0.206 


40 + 0.22 
"•^'^-0.19 


Q1+0.04 
-^-0.03 


1.4 10^5 


23 


A773 


0.217 


7 41 +0-33 


n 2Q+0-05 
^•^^-0.05 


1.5 10^^ 


15 


A 1763 


0.223 


^ ^^A-V\ .34 


^ „ . _i_n Oft 


1.7 10^^ 


12 


A2261 


0.224 




o-4ilg:ii 


2.5 10^^ 


4 


A267 


0.231 


-0.37 


^ „ . _i_n 07 
^•'^^-0.06 


1.2 10^^ 


15 


A2390 


0.231 


^ ^^-l-D 9Q 

8.68lg;2? 


'-'•'^^-0.04 


4.0 10^^ 


12 


RXJ2129. 6+0005 


0.235 


_ _ . _|_n 04 

5.74lg;;^ 


^ ^r^-l-0 0.3 

0.381^3 


2.0 10^^ 


43 


A1835 


0.253 


_ ^ . _i_n 1 n 
7 1 4T^'-'.iu 

' •-'-^-0.11 


'^•'^'^-0.02 


4.9 10^^ 


37 


RXCJ0307.0-2840 


0.253 


_ . _-i-n .38 

-0.35 


^ ^^-1-0 06 


1.6 10^^ 


11 


E1455+2232 


0.258 


^•^^-0.09 


^"^^-0.03 


1.9 10^^ 


33 


RXCJ2337.6+0016 


0.273 


_ _ . _i_n fifi 
' • '^-0.52 


^ _ ^-l-D 07 

n 1 q^^-^' 
^•^^-0.07 


1.7 10^^ 


11 


RXCJ0303.8-7752 


0.274 


91 +0.64 
^•^-^-0.62 


n 2(S+0-08 
'-'•^"-0.07 


1.8 10^^ 


10 


A1758 


0.279 


Q 1 (V^^--^^ 
y--LO_o.43 


n 2Q^ 


1.9 10^^ 


44 


RXCJ0232.2-4420 


0.284 


7 1 Q + 0-31 
' •-^'-*-0.29 


Qn+0.05 
^•'^'^-0.05 


2.6 10^^ 


10 


ZW3146 


0.291 


n r,-, +0.14 

"•^-•--0.10 


n qo+0.02 
"^•'^'^-0.02 


4.1 10^^ 


50 


RXCJ0043. 4-2037 


0.292 


„ „_-|-0 46 

6.95lg:t5 


^ r^r^-l-O 08 

n 2Q^ 

"^•^^-0.07 


1.8 10^^ 


10 


RXCJ0516.7-5430 


0.295 


00+0. 84 
^•'^^-0.74 


n 1Q+0-09 

^•-^^-0.09 


1.7 10^5 


10 


RXJ0658-55 


0.296 




Q 90+0.03 
"^•^'^-0.03 


5.7 10^5 


29 


RXCJ2308.3-0211 


0.297 


' -^^-0.67 


n 40^ 
'^•^^-0.11 


1.3 10^^ 


9 


RXJ2237.0-1516 


0.299 


•^•^"-0.42 


4fi+0-24 
"^•^"-0.20 


3.1 10^4 


19 


RXCJ1131. 9-1955 


0.307 




0.28l°;°e 


2.4 10^^ 


11 


RXCJ0014.3-3022 


0.308 


8.36l°;« 




3.1 10^^ 


13 


MS2137-23 


0.313 


-0.19 




2.1 10^^ 


11 


MS1208.7+3928 


0.340 


r 0^ + 1.04 

^•^^-0.80 


0.79l°|^ 


3.9 10^4 


11 


RXJ0256.5+0006 


0.360 




47+0-18 
-0.17 


9.7 10^4 


11 


RXJ0318. 2-0301 


0.370 


6 07^ 

-0.81 


n 22^ 


8.6 10^4 


16 


RXJ0426. 1+1655 


0.380 


'^•^^-1.73 


0.54t0;33 


8.0 10^4 


10 


RXJ1241. 5+3250 


0.390 


^0+0 74 
6.63t^:^^ 


00+0 15 


7.5 10^4 


17 


A851 


0.406 


0-^^-0.45 


n 22+°-°^ 

'^•^^-0.07 


9.9 10^4 


43 


RXCJ2228+2037 


0.412 


Q QQ + O.SO 


n 90+0.06 

"^•^'^-0.06 


2.6 10^^ 


23 


RXJ1347-1145 


0.451 


11 44+0-26 


PI 97+0.03 
-0.03 


1.1 10^^ 


32 


CL0016+16 


0.541 




29+°-°^ 

^•^^-0.06 


4.0 10^^ 


28 


MS0451. 6-0305 


0.550 


10 02+°-^° 


Sfi+0-08 

'^•'^"-0.07 


4.1 10^^ 


28 


RXJ1120.1+4318 


0.600 


6 09+°-^^ 

"•"^^-0.69 


'-'•^^-0.16 


1.3 10^^ 


18 


MS1137.5+6625 


0.782 




S5+0-37 


1.6 10^^ 


18 


MS1054.4-0321 


0.823 


^•^^-1.03 


"^•^^-0.13 


3.0 10^^ 


25 


WARPJ0152.7-1357 


0.837 


7 QQ+0.73 


"^•^^-0.11 


2.0 10^^ 


48 


CLJ1226. 9+3332 


0.890 


10.6912-82 


15+0-08 
^•-'-^-0.08 


4.2 10^^ 


69 



Note. — Redshifts, temperatures, metal abundances (wrt. Solar), luminosities and exposure times for 
the sample. Redshifts are taken from the NED listing. Temperatures and metal abundances were derived 
using spectral fits fSection lS.lj) . Luminosities within a 1 Mpc aperture were derived using the SPI modeling 
(Section [3:21. 
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3.1. Standard analysis 

All clusters are first analyzed using "standard" spec- 
tral analysis. This was conducted by extracting X-ray 
counts from a circular region centered on the peak of X- 
ray emission. The extraction radius was determined by 
estimating the radius at which a circle encompasses 90% 
of the background subtracted surface intensity. 

The background was estimated using the surface in- 
tensity at the edge of the field. Background spectra 
were extracted in regions outside the source extraction 
region. Point sources were detected using SAS routine 
emldetect using only sources measured with a likelihood 
above 100. 

The e x tracte d spectra were fitted using XSPEC 
(lArnaudl I1996D s oftwa r e, emp l oying a MEKAL 
(I Me we et al .1 I1985L 119861 : iKaastral 1 19921 : iLiedahl et al.l 
1 19951 ) thermal plasma model with Solar abundances 
absorbed by a WABS ([Morrison fc McCammonI Il983[ ) 
model, which we allow to be fitted as a free parameter. 
In the fit, the redshift was fixed to the known optical 
value as listed in the NED database entry. 

The results of these "standard" spectral fits including 
plasma temperature and metal abundances w.r.t. Solar 
are listed in Table [H The bolometric luminosities shown 
in the table are derived using the SPI analysis below. 
The redshift as given in NED is also shown as well as 
the average effective exposure time after filtering. All 
observations had usable data for all three EPIC detectors 
with the exception of A665, A1413, A2261, A2597 and 
A3921 where only the two MOS detectors were available. 

3.2. SPI based MCMC analysis 

The cluster event files are modeled us- 
ing Smoothed Particle Infe rence (SPI) 

([Peterson. Marshall, fc Anderssonl I2007D with a Monte 
Carlo model of the XMM-Newton EPIC Camera 
(|Andersson et all l2007f ). Within the SPI analysis, 
clusters are modeled as conglomerations of 2D spatial 
Gaussians with individual spectral models. In this work 
we use the MEKAL model to describe the thermal 
plasma. Model photons are simulated and propagated 
through the detector model adding background Monte 
Carlo events representing internal fluorescent lines, 
electronic noise and soft proton signals. Data and model 
photons are binned in three dimensional adaptive bins 
and compared via a two-sample likelihood function. 
All parameters, spatial and spectral, are iterated in a 
Markov Chain with an adaptive step length. All model 
samples in the converged part of the chain are used to 
reconstruct the cluster. 

3.2.1. SPI setup 

The number of SPI particles are determined by the 
number of photons in the data so that, on average, 
N^/Np = 400, where Np is the number of particles and 
is the number of X-ray events in the data. This 
number i s chosen based on the Bayesian evidence calcu- 
lation in iPeterson. Marshall, fc Anderssonl (|2007l ) where 
we find the evidence reaching a plateau near 400 SPI 
particles for a 155000 photon observation. We choose to 
scale the number of particles with the number of photons 
since that is what ultimately determines the complexity 
of the data and we want a corresponding complexity of 



the model. We constrain the minimum number particles 
to be 100 and the maximum to be 1000. 

The oversimulate factor, the factor that determines the 
number of model photons, is set so that there are 10 
times as many photons in the model compared to the 
data with a maximum of 4 million model photons. This 
number is motivated by the drastic improvement of the 
opti mization when using a factor of 10 or ab ove as shown 
in Peters on. Marshall, fc Anderssonl (|2007l ). The upper 
limit is set to minimize computing time for very well 
exposed clusters. 

The three dimensional (x,y,pulseheight) adaptive bin- 
ning grid is created so that every bin with more than 
20 photons is divided in 2. The spectral dimension is 
divided 10 times more often than the spatial dimensions 
on average in order to achieve appropriate spectral reso- 
lution. 

3.2.2. Model setup 

The setup of the instrumental background model, con- 
sisting of electronic noise, soft proton detections and 
fluorescent emiss i on lin es, is analogous to the setup in 
lAndersson et al.l (|2007l ). The fraction of photons going 
to the background model is variable from to 1, as are 
the relative normalizations of the included components. 

We model the soft X-ray background originating from 
our Galaxy using a uniform emission component con- 
sisting of a thermal plasma spectral model. This emis- 
sion consists of local (and thus weakly absorbed) compo- 
nent plus an absorbed, more distant contribution from 
the Galactic halo; it is adequately described as sev- 
eral thermal component s in the 0.05-0.5 keV range (e.g. 
iKuntz fc Snowdei3l2000l ). Motivated by our analysis of 
blank sky data flies, we flnd that unabsorbed MEKAL 
model at temperature of 0.16 keV, with metal abundance 
of 0.3^0 at z = describes the data well, and this is 
a model we adopt for the Galactic background: given 
the limited EPIC bandpass, our approximation is suffl- 
cient. We keep in mind that the spectral accuracy of the 
method is limited in the regions of the clusters where a 
large fraction of the photons come from the X-ray back- 
ground, such as at large radii. While there the cluster 
flux is low, the adaptive binning grid will allow only gross 
spectral features to be detectable. 

Similarly the Cosmic X-ray Background (CXB) - pre- 
sumably due to superposition of unresolved AGN - is 
modeled using a powerlaw with photon index F = 1.47 
and absorption fl xed at the galactic v alue at the coordi- 
nates as given bv lDickev fc LockmanI ([1990). 

The cluster model consists of spatial Gaussians de- 
scribed by an x and y position and a Gaussian a. Each 
Gaussian is assigned a spectral MEKAL model with 
WABS absorption. The allowed ranges for the absorb- 
ing equivalent hydrogen column, uh^ the plasma tem- 
per at ure^_^^^_themetall^ w.r.t. the Solar values of 
Ande rs~fc Greves"s3 (| 19891 ), Z, and redshift z are shown 
in Table [3l Absorpti on is assumed to b e with in 20% 
of the riH.gai values of iDickev fc LockmanI (|1990[ ) at the 
cluster coordinates. 

3.2.3. Convergence 

In IPeterson, Marshall, fc Anderssonl (|2007l ) a criterion 
for convergence of the Markov chain is applied and it 
is found to converge within 200 iterations. Here, the 
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TABLE 3 

Cluster model parameter ranges used 
in the spi analysis 



PcircirQGtGr 


min value 


max value 


Spectral model 






riH 


0.8 nH,gal 


1-2 nH,gal 


T [keV] 


0.5 


15.0 


Z wrt. Solar 





2 


z 


ZNED fixed 




Spatial model 






X 


-12^ 


+12' 


y 


-12' 


+12' 


In a{") 
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Note. — The parameter ranges used in the 
SPI modelin g, uh is variable wit hin 2 % of 
the values of |Dicke\^ Lockma^ (llQQOT) . The 
X and y positions are variable within 12' of the 
XMM nominal pointing. 



more conservative limit of 750 iterations is used as the 
point of convergence. All the results derived from the 
model samples are from the iterations from 750 to 2000. 
This range is chosen based on slight deviations in cluster 
properties when derived from iteration 200 and onward. 

3.2.4. Simulated clusters 

In order to study the systematic effects induced by the 
X-ray mirror PSF and limited statistics on faint sources, 
we simulate a massive, spherically symmetric, weak cool- 
ing core cluster at z = 0.1, 0.2, 0.4 and 0.8. The simu- 
lated cluster model consists of two superimposed beta 
models with core radii of 74 and 370 kpc, both with 
/3 = 0.7. The smaller, core component is assigned a 
spectral model with kT = 4 keV whereas the ambient 
component has kT = 9 keV. The normalization ratio of 
the cold to the hot component in terms of emission mea- 
sure, J nefipdV^ is 11/9 and the overall normalization is 
set so that the total bolometric luminosity of the cluster 
is LboI = 1-7 10^^ erg s~^, a typical luminosity in our 
sample, present at all redshifts 0.1 < z < 0.8. The model 
is based on observations of the cool core clusters A2029 
and A2241. Bolometric luminosities for the cluster sam- 
ple are estimated using the 0.01 - 100 keV energy band 
as described in Section [4?T1 

These simulated clusters all have the same background 
level (5 10~^ s~^ assumed to have an effec- 

tive exposure time of 20 ks in all 3 EPIC detectors. The 
assumed luminosity results in a total number of 850058, 
209294, 46358 and 9353 cluster photons respectively for 
the z = 0.1, 0.2, 0.4 and 0.8 clusters in the adopted en- 
ergy range (MOS:0.3-10 keV, pn:l.l-10 keV). As a refer- 
ence, the same clusters are also simulated without back- 
ground and reconstructed separately. 

The apparent evolution with z in the derived properties 
for these identical clusters will be used to categorize the 
clusters in the sample. Figure [2] shows the original cluster 
model before propagated through the instrument model 
visualized with a luminosity and temperature map within 
the 1 Mpc radius. Luminosity per unit solid angle in 
this Figure is given in units of 10^^ erg s~^ {")~'^ and 
temperature in keV. 

In order to follow the evolution of substructure in clus- 
ters we also simulate two sets of clusters with irregular 
morphology, one using a two "subcluster" model and the 



other, a three "subcluster" model. The two-subcluster 
model consists of two beta models, at 4 and 9 keV respec- 
tively, each with Tc = 74 kpc and /3 = 0.7 separated by 
300 kpc. The total luminosity of the clusters is the same 
as for the original cool core cluster, LboI = 1-7 10^^ erg 
s~^. The cluster model is shown in Figure [3l The three- 
subcluster model consists of three beta models, at 4, 9 
and 9 keV respectively, with core radii of Tc = 74 kpc 
and (3 = 0.7. The clusters are separated by a triangle 
with sides of 300, 200 and 200 kpc. The total luminosity 
of the three clusters is LboI = 1-7 10^^ erg s~^. This 
model is shown in Figured! 

4. POST-PROCESSING 

Here, we describe the raw output of the SPI analysis 
- temperature and luminosity maps of the clusters con- 
sidered here - and discuss the estimation of uncertain- 
ties. We also outline the methods that will be applied 
in the subsequent sections towards determination of the 
luminosity and temperature structure of clusters in our 
sample. First, we consider a luminosity and tempera- 
ture contrast analysis with the primary goal to identify 
clusters containing cool co res. The luminosity con trast 
analysis is similar to that of IVikhlinin et all (|2006bl ) . We 
also co nsider t he ^T ower Ratios" method, suggested by 
iBuote fc Tsai ( 1995f ) to quantify the cluster substructur- 
ing. Finally, we propose a new method, the "Tempera- 
ture two-point correlation" which is specifically designed 
to quantify the temperature structure, and is enabled by 
the SPI analysis. 

4.1. Temperature and luminosity maps 

We create 2x2 Mpc cluster luminosity and emission 
weighted temperature maps centered on the peak of clus- 
ter emission with 10 kpc bins. The me thod of creating 
median parameter maps is described in Andersson et al\ 
(|2007l ). Maps are created for each sample in the chain 
from iteration 750 to 2000 and are averaged by taking 
the median in each spatial point over the whole sample. 

Point sources are removed by filtering out all cluster 
particles within a radius of — 16.1 + logio(^2) 15 '^ from 
the point source as detected by emldetect, where L2 is 
the likelihood of detection. We find that this method suc- 
cessfully removes any point source contamination with- 
out removing the cluster flux in the region of the point 
source. This method is effective because particles that 
represent the point source emission are generally smaller 
in size whereas the cluster particles are larger and fill the 
region where the point source was removed. An example 
of such point source removal is shown in Figure [5] which 
shows maps of luminosity (top) and temperature (bot- 
tom) for Abell 3888 before (left) and after (right) point 
source removal. 

The reconstruction of the four simulated cool core clus- 
ters within 2x2 Mpc is shown in Figure [6] where lu- 
minosity per unit solid angle is given in units of 10^^ erg 
s~^ ('0~^ temperature in units of keV. Here, the ef- 
fects of the XMM point spread function (PSF) as well 
as the loss of photons with redshift is clearly seen as a 
distortion and flattening of the profile at high z. 

The maps for a subsample of clusters are shown in 
Figures [16] through [32] including the named clusters in 
Figure [H We comment briefly on these maps in Section 
15.41 and compare them to previously published results. 
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Fig. 2. — 2x2 Mpc images of the bolometric luminosity in units of 10^^ erg s ^ {") ^ (left panel) and temperature in keV (right panel) 
of the original simulated cool core cluster model. 




Fig. 4. — 2x2 Mpc images of the bolometric luminosity in units of 10"^^ erg s ^ ^ (left panel) and temperature in keV (right panel) 
of the simulated three-subcluster model. 



Note that the temperatures in Table [T] do not necessarily 
agree with the scale shown in the figures. This is because 
these values are derived using two different methods. The 
tabulated values are determined using standard analysis 
(Section 13. ip requiring a single value of gas temperature 
for the whole cluster, whereas the parameter maps are 
created using a range of temperatures derived from the 
SPI runs (Section [3T2|) . The source of the disagreement 



has to do with the fact that SPI bases the spectral mod- 
eling on superpositions of thermal spectra with many 
different temperatures. The fact that the similarities 
between these spectra increase with higher temperature 
implies that the average temperature value of these also 
increases. It also becomes dependent on the prior range 
of temperatures (here: 0.5 - 15 keV). Whenever referring 
to the cluster average temperature, the tabulated values 
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Fig. 5. — 10^ x 10^ field showing luminosity (top) and temperature (bottom) maps of Abell 3888 before (left) and after (right) point 
source removal. The color scale in the luminosity map is set so that white corresponds to the maximum cluster fiux. The point source is 
100 times brighter than this level. The scale in the temperature map ranges from 2 to 10 keV. 
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Fig. 6. — Luminosity (top) and temperature (bottom) maps for the reconstructed simulated cool core clusters, with background, at 
redshift 0.1, 0.2, 0.4 and 0.8. Luminosity per unit solid angle is in units of 10^^ erg s"-*^ and temperatures in keV. 



are used. A bolometric correction is applied to the clus- dance and emission measure of the particle are used to 
ter particles individually, where the temperature, abun- calculate a bolometric luminosity in the 0.01 eV to 100 
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keV range using the MEKAL spectral model. This is the 
luminosity that is listed in Table [H 

4.2. Estimation of uncertainties 

In all cases, except for the temperature and elemental 
abundance values used in Table (TJ uncertainties are esti- 
mated by identifying 5 independent sets of samples in the 
Markov chain posterior. The samples are taken sequen- 
tially in ranges containing 250 iterations each, starting 
at iteration 750. The luminosity and temperature maps 
and the quantities derived from these are then evaluated 
separately and the uncertainty on these quantities is de- 
termined from the rms scatter within these 5 sets. This is 
a valid approximation of statistical uncertainty because 
each subset can be seen as an independent reconstruction 
of the cluster. The motions of particles in the parameter 
space are large and within a few iterations there is no 
memory about earlier states. This is largely due to the 
high dimensionality of the overall parameter space and 
the lack of local likelihood maxima. Of course, it would 
be more accurate to use more than 5 sets for this calcu- 
lation but at least 250 iterations are needed in order to 
describe a cluster model and there is a limitation in CPU 
time to produce more iterations. 

4.3. Luminosity and temperature contrast methods and 
identification of cooling clusters 

A common definition of a "cooling fiow" cluster regards 
the central cooling time, tcooi = T/d{ln{T)), being much 
less than a "Hubble time", tcooi <<tH- The calculation 
of the cooling time requires high resolution spatially re- 
solved spectroscopy which is not available for most of the 
clusters in our sample. Instead we use the information 
about the steepness of the luminosity profile as well as 
the gradient of temperature in the cluster core to empir- 
ically select the clusters with cool cores. 

The luminosity contrast is calculated as 



Cl 



L{r< 0.02r5oo) 



L(0.08r5oo <r< 0.12r5oo) 



(1) 



where L is luminosity per solid angle. We use the lumi- 
nosity based definition of rsoo; 



rsoo = 909 E{z) 



^bol 



lO^^erg s- 



0.172 



kpc. (2) 



where 



E{z) = H{z)/Ho = 
+ z)2 {I^Qm z)-z {2^ z) Qa. (3) 

This value of r5oo is de rived by combining the M-T 
relat ion fromlArnaud et al . (2005) and the L-T relation 
from lArnaud fc EvrardI (| r999). The fractions of rsoo are 
used to account for the difference in size for clusters of 
various luminosities. 

Similarly we define the temperature contrast as 



Cj 



T (r < 0.02r5oo) 



T (0.45r5oo <r< O.SSrgoo) 



(4) 



where we choose to estimate the gradient farther from 
the core (0.5r5oo) than for the luminosity contrast above 



because it is less sensitive to smearing by redshift depen- 
dent effects. We find that the radius of 0.5r5oo is where 
we could get the most distinguishing power. We discuss 
the use of Cl and Ct towards identifying cooling clusters 
in Section 

4.4. Power Ratios method 

One of the successful approaches towards assessing 
the amount of substructure i n clu sters is the power 
ratios method (iBuote fc Tsail Il995l ) . It was used in 
iJeltema et all (j2005[ ) for a sample of Chandra-observed 
clusters to establish that clusters are more dynamically 
active at high z. This method is based on the multipole 
expansion of the surface brightness (/)) around the 
cluster centroid where the multipole moments and bm 



am{R)= ^{R,(I)){R')'^ COS m(l)'d^x' (5) 

Jr'<r 

bm{R)= [ T.{R,(})){R')'^smm(})'d^x' (6) 
Jr'<r 

so that the powers in the multipole m can be written 



Po = (aoln(i?))' 



1 



2m?R^'^ 



(7) 
(8) 



Here, we define the location of the cluster centroid by 
the requirement that Pi should be zero at this location. 
The above method is applied to the luminosity maps de- 
scribed in the previous section and the ratios P2 / Pq and 
P3/P0 are calculated within a radius of 500 kpc. We 
expect that P2/P0 will be larger for elongated clusters 
whereas P3/P0 will be large for clusters with much sub- 
structure in the luminosity map. 

4.5. Temperature two-point correlation 

We suggest that the variation of temperature inside the 
clusters can be estimated and characterized by taking a 
(non-standard) two point correlation of the temperature 
difference, weighted by products of luminosity so as to 
enhance temperature differences among the regions with 
the highest brightness. This is accomplished by taking 
the following sum over all 10 kpc x 10 kpc pixels in the 
generated maps: 



A{rk) 



(9) 



where Li is the luminosity of pixel Ti is the tem- 
perature and Lt = y^Pi Lj. i and j are such 
that both {xi^yi) and (xj^yj) cover all points within 
500 kpc of the centroid calculated in Section 14.41 A 
is binned based on the distance between i and j as 
10k kpc< rk < I0{k ^ 1) kpc, k = 0...99, where 

r/c = ^y {xi — XjY + ^Vi — V3Y is pixel distance and 
(xi^yi) is the spatial position of pixel i. 

This statistic is designed to quantify the amount of 
distortion in the intra cluster medium, specifically re- 
garding temperature features that can result from clus- 
ter mergers. Basically, it gives a "power spectrum" of 
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strong temperature gradients over different scales, where 
the "strength" of the gradient is determined by ^JLi Lj. 
In integrating this quantity over smah or large distances 
it should be possible to distinguish small scale features, 
such as cooling cores, from larger scale disturbances. We 
discuss the application of this method in Section 15.3.11 

5. RESULTS 
5.1. Identification of cooling core clusters 

The first step is to identify clusters with cooling cores. 
We use the approach described in Section 14. 3[ and plot 
the luminosity contrast vs redshift in Figure [3 (left 
panel). Here, we also show the names for some well 
known clusters. Low Lx clusters (< 4 10^^ erg s~^) are 
shown as red circles, intermediate Lx clusters (4 10^^ erg 
s"-*^ < Lx < 2 10^^ erg s~^) as green stars and high Lx 
clusters (> 2 10^^ erg s~^) as blue squares. Likewise, 
the temperature contrast vs redshift is shown in Figure 
[71 (right panel). 

In Figure [71 we also show the calculated values for our 
simulated cool core cluster at redshift z = 0.1, 0.2, 0.4 
and 0.8. It is easily seen that the effects of PSF smooth- 
ing and loss of photons due to increased distance severely 
affects both quantities. We show both the simulated clus- 
ters with (dark gray) and without (light gray) simulated 
background, totaling 8 simulated clusters. 

Cl is estimated at 0.1 rsoo (~ 10'^ at z = 0.8) and it be- 
comes increasingly difficult to estimate at high redshifts. 
This can be seen in the sharp drop in Cl in the trend 
of the simulated clusters above z = OA. At 2; = 0.8, the 
detection of cooling cores using Cl is no longer sensitive 
and thus we decide to rely on the Ct statistic alone for 
identification purposes above z = 0.6. We identify those 
clusters with a Cl above a line interpolated between the 
points set by the simulated clusters as cool core clus- 
ters. This cut corresponds to a linear interpolation of 
the points Cl = (7,7,5,2.2) at z = (0.1,0.2,0.4,0.6). 
This comfortably separates known cool core clusters (e.g. 
A2029) from disturbed and intermediate core clusters 
(e.g. A1689). 

For Ct, the simulated clusters can be seen to have a 
trend that is more linear and not as dramatic as the trend 
in Cl- This is largely due to the larger radius (0.5 rsoo) 
used when calculating Ct- Using Ct to identify cool 
core clusters we make a cut corresponding to a straight 
line approximately following the trend of the simulated 
clusters so that clusters below the line, Ct < 0.85 + 
0.15 logz, are selected. These cuts select well known 
cooling clusters like A1068, A1835 and A2204 but fail to 
select clusters that are known not to have pronounced 
cool cores such as A1689 and A1413. 

In subsequent plots we show the values for the recon- 
structions of the simulated clusters as filled circles and 
use this to measure how accurately a weak cooling core 
cluster at different redshifts can be resolved. The cuts 
are shown in Figure [3 In Figure [51 we also show Ct plot- 
ted against Cl, clearly showing the separation of cool- 
ing clusters like A1835 and A2204 (G, J - shown bot- 
tom right), nearly isothermal clusters such as A1413 and 
A1689 (B, K - shown center) and disturbed clusters like 
A2218 and A520 (H, E - shown top left). 

Out of 101 clusters, 31 are identified as cooling clusters. 



5.2. Dynamical activity 
5.2.1. Power ratios 

In assessing the dynamical activity in the clusters in 
our sample we perform a multipole expansion of our lu- 
minosity maps as described in Section 14. 4[ The power 
ratios P2 / Pq and P3 / Pq are calculated to a radius of 500 
kpc in order to quantify the amount of substructure in 
the sample. These are listed in Table [H 

Here we divide the sample in 3 redshift bins; a local z 
bin {z < 0.1), a low-z bin {0.1 < z < 0.3) and a high-z 
bin {z > 0.3). The results of the multipole expansion are 
shown in Table [6l where we show the average values of 
P2/P0 and P3/P0 for the 3 samples. We also show the 
result obtained when using only the cooling core clusters 
identified in previous sections. These are not excluded 
from the other samples. 

The values of P2/P0 and P3/P0 are plotted against 
z in Figures [9l and [TOl respectively where the simulated 
cool core clusters are shown as filled circles connected 
by black lines. In Figure [H we also show the simulated 
two-component clusters as purple triangles connected 
by black lines and in Figure [TOl we show the simulated 
three-component clusters as yellow triangles connected 
by black lines. 

The simulated cool core clusters represent the min- 
imum amount of substructure that can be measured 
as they are simulated with perfect circular symmetry. 
This is true for clusters of similar luminosity and expo- 
sure. However, the P3/P0 values in Table [HI are higher 
than those of the simulated cool core clusters and above 
2: = 0.3 most clusters have P3/P0 values similar to those 
of the simulated three-component clusters. This is most 
clearly seen around z = 0.4: where the the datapoints 
cluster around the value of the three-component simu- 
lation. Towards z = 0.8 the simulated irregular clus- 
ter shows a decreasing value with redshift due to the 
smoothing effect introduced by PSF blurring and lack 
of photons. The datapoints around z = 0.6 and z = O.S 
show large Ps/Pq values around 10~^, larger than for the 
simulated clusters. 

To investigate further the proposed evolution of P3 / Pq 
we restrict the analysis to clusters below z = 0.6 due 
to the apparent breakdown in the analysis as seen in 
the simulated clusters at z = 0.8. We also exclude the 
low-luminosity clusters (P50Z < 4 10^^ erg s~^) from the 
sample since these are not present at all z. We bin the 
P3/P0 values for the remaining clusters in 5 bins of red- 
shift {z = 0,0.1,0.2,0.3,0.4,0.6) and fit these datapoints to 
a straight line in log z-log P3 /Pq space. The binned data 
is shown along with the best fit line in Figure [Til In 
order to distinguish between real and apparent evolution 
we interpolate the slope linearly for values of log P3 / Pq 
from the simulated cool core clusters (gray circles) to 
the simulated clusters with irregular morphology (yellow 
stars). These interpolated lines are shown as dashed lines 
in Figure [m where the red dashed line has the lowest 
when compared to the data. Compared to the best line- 
fit (solid line) the dashed line models have Ax^ = 2.6, 1.1 
and 3.3 respectively (from bottom to top). We conclude 
that the evolution in P3/P0 in addition to the redshift- 
dependent bias is only slightly more than 1-a significant. 

To compare our results with pre vious work we l i st our 
results for the clusters also found in lJeltema et al.l (|2005l ) 
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Fig. 7. — Luminosity contrast (left) and temperature contrast (right) as a function of redshift. The different symbols represent luminosity 
(red circles > 10^^ erg s"-"^, green stars > 4 10^^ erg s"-*^, and blue squares > 2 10^^ erg s"-*^). The cuts on Ci, and Ct are shown as 
solid lines. The cooling clusters are above and below the line respectively. The simulated cool core clusters without (with) background are 
shown as light (dark) gray filled circles, occasionally overlapping. A few selected clusters are denoted using letters from A to O as follows: 
(A) RXJ0658-55, (B) A1413, (C) A1750, (D) A1775, (E) A520, (F) A2029, (G) A1835, (H) A2218, (I) MS1054.4-0321, (J) A2204, (K) 
A1689, (L) RXJ1347-1145, (M) ZW3146, (N) MS2137-23, (O) A1068. 




Fig. 8. — Temperature contrast plotted versus luminosity contrast. Here the different symbols represent redshift (red circles 0.069 < z < 
0.1, green stars 0.1 < z < 0.3, and blue squares z > 0.3). Simulated cool core clusters are shown as filled gray circles. Letters from A to O 
denote a few selected clusters (see Figure (T)). 



along with their results in Table [3 These values were de- 
rived using Chandra data only. We find our values to be 
consistently about a factor 2 higher with a few excep- 
tions but that they follow the same trend of low to high. 
This discrepancy could be related to the limited number 
of photons in the Chandra data causing structures to ap- 
pear smoother but could also be due to artifacts in our 
modeling described further in following sections. 

5.3. Luminosity- Temperature relation 

Using our definition of a cooling cluster, we treat these 
clusters separately in the subsequent analysis. The re- 
maining clusters are divided into 2 redshift bins; 0.069 < 
z < 0.2 and z > 0.2. This results in 38 low redshift 



clusters with z = 0.11 and 32 high redshift clusters with 
z = 0.39. 

The X-ray luminosity-temperature relations for these 
3 samples are shown in Figure [121 We have used the 
bolometric luminosity as calculated within a 1 Mpc ra- 
dius along with the temperatures given in Table [TJ The 
low redshift sample (red circles), the high redshift sam- 
ple (green stars) and the cooling cluster sample (blue 
squares) are shown along with their best fit Lx — T re- 
lations (solid lines). The best fit relations from the n on- 
cooling and cooling samples of lAllen fc Fabiar] (|l998l ). at 
z = 0, are shown as dashed lines for comparison. 

Our best fit values of Lq and a^r are shown in Ta- 
ble [HI where the bolometric X-ray luminosity is Lx = 
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TABLE 4 
Power ratios 



Name 


P2/Po(xlO-^) 


P3/Po(xlO ^) 


CTZA T1fi4^ 4_7SS4 


^ . 1 Z/ _i_ . c/tj 


n 41 _L n 4q 


A1837 


112 zb 7 


9 q7 _|_ n qi 
Zi.O I ZIZ u.o± 


A3112 


on c: 1 1 


n fjQc: _L n nqi 


A1775 




q Of? _|_ fi 
O.OU ZIZ u.uo 


A399 




q p:?: _|_ n qi 

0.00 ZIZ U.O-L 


A1589 


298 zb 17 


n QQ _L n cr-i 
vj.oo _i_ yj-Kji. 




208 zb 6 


U.UOc/ ZIZ U.UOZi 


A401 


104 zb 3 


n q9 _L n 07 


A2670 


c/U ZIZ U. c/ ± 


4 07 -1- 4*^ 

^.U ( ZIZ u.^o 


A2029 


Q _|_ 9 1 
tJO. c/ ZIZ J. 


n nc:i _|_ n 099 




91 _L 
Z/ ± . u _i_ . \J 


Q q7 _L 9 79 

. c/ 1 _i_ Z/ . 1 Z/ 




78.0 zb 7.8 


27.1 ± 1.8 






n c:o _L n 49 


"RYn 70^91 «-|-ni19 


106 zb 16 


n Q7 _|_ n CO 
u.o ( ZIZ u.oo 




45.0 zb 11.7 


n c^q _|_ n f^q 
U.Oc/ ZIZ u.uo 


A1650 


07 q _L 9 q 


n 94 _L n Of? 

U.Zf^ _l_ u.uu 


A1651 


00. U ZIZ 0. ± 


n 1 9 -u OS 

u. ±Zi ZIZ u.uo 


A2597 


9(^ _L n 

Z/\J . _I_ L/. \J 


U.UOO _l_ U.UZ/U 


A1750 


qq q _l_ c; 7 

Oc/. c/ ZIZ . 1 


1 04 _|_ fl 7'^ 

± . 0^ ZIZ U. ( 


A478 


4/^ 4 _|_ n c; 

^U.^ ZIZ (J.O 


n nc:4 _|_ n 004 
U.UO^ ZIZ u.uu^ 


A278 


16.2 zb 2.1 


1 9q _L n c:7 


A2142 


247 zb 9 


u. ±0 ZIZ u.uo 


A3921 


990 -1- 1 n 

Zij^iO ZIZ -LU 


9 9 _|_ fl 48 
^.Zi ZIZ u.^o 


A13 


00 9 _|_ 1 n 1 
00. Zi ZIZ ±u. ± 


fl 1 q -U 1 9 

U. ±c/ ZIZ \J . LZi 


A3911 


344 zb 11 


9 4c _|_ fl CI 
Z.^O ZIZ U.O-L 


T^xr! 79S1 Q fi-7S1 S 


128 zb 18 


fl 99 _L fl nq 

U.ZfZf _l_ U.Uc/ 




19 c: _|_ q 

±^.0 ZIZ c/.U 


99 q _|_ f: c 

ZiZi.c) ZIZ U.o 


A3827 


17.0 zb 1.7 


fl qfi -U 1 ^ 

U. c/U ZIZ U. -LO 


7^Xn7n911 4 4017 


18.3 zb 1.2 


1 40 _|_ fl CC 
L .^0 ZIZ U. 00 


A2241 


7 ni _|_ 1 qq 

( . U J- ZIZ J- . Oc/ 


fl 40 _|_ fl qf? 

U.^O ZIZ U.OU 




36.0 zb 1.7 


fl ni Q _L fl ni 
u.u±o _l_ u.u± 


pxn 7nfi4^ 4-^41 


126 zb 6 


911 _|_ fl f?7 

Zi. L L ZIZ U.U ( 




1 q q _|_ 9 q 

J-O.O ZIZ Zi.O 


fl 1 q _|_ fl 14 

U. -Lc/ ZIZ U. -L^ 


A1302 


9Q q _|_ q n 

AiO.a ZIZ O.U 


fl 9f? _|_ fl 9q 
U.^U ZIZ u.^o 


7^Xn7nfi1fi R-474R 
J uuxu.o-^ ( 


112 zb 18 


3.0 ± 1.7 


7^Xr!7914Q 1-^041 


9q q _L q 7 


fl fi90 _L fl ni c 
U.UZ/O _i_ u.u±o 


7^xn7if^ifi "^-i-nnnf^ 


09 q _|_ 4 n 

OZi.a ZIZ ^.U 


fl f^q _|_ fl 41 
u.uo ZIZ u.^± 


RXP 71 1 41 4 1 91 fi 
ixyvv^ J J. J. ^ J. . J. z J. u 


19.4 zb 1.5 


fl 1 9 _|_ fl 1 1 

U. J-Z ZIZ U. J- -L 


7^Xn 70090 7-9'=;49 


118 zb 6 


1 qf? _|_ fl qi 

L . OU ZIZ U. ± 


7^Xn71044 "=^-0704 


47.0 zb 3.0 


fl 77 _|_ fl nq 

U. ( ( ZIZ U.Uc/ 


RXCJ0145.0-5300 


342 ± 29 


2.63 ±0.93 


A1068 


63.0 ±4.4 


1.21 ±0.32 


RXJ1416.4+2315 


180 ± 20 


3.37 ±3.16 


RXCJ0605.8-3518 


47.0 ±3.0 


0.52 ±0.17 


A1413 


192 ±3 


0.19 ±0.1 


RXCJ2048. 1-1750 


34.4 ±4.4 


3.92 ± 0.34 


A3888 


117 ±7 


6.33 ±0.8 


RXCJ2234.5-3744 


119 ±4 


4.06 ± 0.49 


A2034 


84.5 ±8.4 


0.41 ±0.21 


A2204 


3.81 ±0.33 


0.077 ± 0.045 


RXCJ0958.3-1103 


79.0 ±9.1 


0.14 ±0.12 



Note. — Results of the Power Ratio analysis (Section 
14. 4|) of the clusters. 

LQ{kT/6 keV)"^^. We have added a systematic scatter 
in log(T) in order to achieve x^/dof. = 1. 

It is obvious in this plot that it is of great importance 
to exclude cooling core clusters since these clusters have 
much higher Lx for a given T compared to non-cooling 
clusters. It is of particular importance to select these 
clusters using a method unbiased in redshift. A selection 
bias at high redshift leading to the failure to identify 
cooling clusters properly could easily be interpreted as 
evolution. We have based our selection of cooling clus- 
ters on simulated identical clusters at various redshift s, 
propagating their photons through our detector model 
to best account for any distance dependent systematic 
effect. 



Fitting the data to a generalized Lx-T relation with a 
z dependence; 

weak evolution is found with the best fit values for the 
non-cooling clusters shown in Table [9] along with fits for 
high luminosity {L^oi > lO^^erg s~^) non-cooling clus- 
ters, high P3/P0 clusters {P3/P0 > 10"'^) and for the 
cooling clusters separately. 

The results for the non-cooling clusters show weak evo- 
lution whereas the cooling cluster sample is consistent 
with no evolution. There is however large intrinsic scat- 
ter within the cooling sample. Interestingly, the sample 
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TABLE 5 
Power ratios 



Name Ps/^oCx lO"'^) P3/^o(x lO"'^) 



A868 

RXCJ2014.8-2430 
A2104 

RXCJ0547.6-3152 

A2218 

A1914 

A665 

A1689 

A383 

A520 

A2163 

A209 

A963 

A773 

A1763 

A2261 

A267 

A2390 

RXJ2129.6+0005 
A1835 

RXCJ0307.0-2840 

E1455+2232 

RXCJ2337.6+0016 

RXCJ0303.8-7752 

A1758 

RXCJ0232.2-4420 
ZW3146 

RXCJ0043.4-2037 

RXCJ0516.7-5430 

RXJ0658-55 

RXCJ2308.3-0211 

RXJ2237.0-1516 

RXCJ1131. 9-1955 

RXCJ0014.3-3022 

MS2137-23 

MS1208.7+3928 

RXJ0256.5+0006 

RXJ0318.2-0301 

RXJ0426.1+1655 

RXJ1241. 5+3250 

A851 

RXCJ2228+2037 

RXJ1347-1145 

CL0016+16 

MS0451. 6-0305 

RXJ1120.1+4318 

MS1137.5+6625 

MS1054.4-0321 

WARPJ0152.7-1357 

CLJ1226. 9+3332 



105 + 6 
15.8 + 1.2 

68.0 + 16.7 
29.4 + 3.1 
72.8 + 5.6 
32.4 + 1.0 
44.8 + 3.5 
27.6 + 2.0 

I. 63 + 0.26 

65.3 + 11.2 
45.6 + 5.9 

98.8 + 5.5 
8.14 + 1.32 
83.2 + 11.2 
218 + 8 

12.9 + 11.4 

95.8 + 9.2 
149 + 7 

66.1 + 1.4 

10.2 + 0.4 
14.6 + 3.8 
16.1 + 0.9 
282 + 29 

37.4 + 13.7 
502 + 13 
50.6 + 4.8 

II. 6 + 1.3 

58.5 + 12.9 
296 + 54 
116 + 3 

23.3 + 4.2 
198 + 41 
142 + 18 
5.54 + 4.01 
1.47 + 0.96 
322 + 293 
6.12 + 2.42 
39.7 + 20.2 

35.9 + 37.5 
27.4+17.9 

304 + 30 
144 + 27 
32.6 + 1.0 
114 + 17 
82.3 + 5.1 
137 + 49 
11.5 + 11.9 
220 + 73 
2880 + 260 
11.1 + 3.2 



0.92 + 0.98 
0.37 + 0.13 
2.83 + 2.54 
0.33 + 0.14 
1.14 + 0.36 
2.45 + 0.26 
7.75 + 0.77 
0.64 + 0.12 
0.62 + 0.17 
3.87 + 0.46 
10.7 + 1.7 
1.18 + 0.83 
1.34 + 0.67 
0.82 + 0.39 
1.36 + 1.3 
2.79 + 1.7 
0.92 + 0.64 
2.97 + 0.86 
0.21 + 0.03 
0.36 + 0.09 
2.23 + 0.51 
0.12 + 0.08 
1.16 + 0.42 
2.01 + 0.87 
0.78 + 0.32 
2.5 + 0.62 
0.37 + 0.08 
2.31 + 2.39 

8.28 + 5.55 
8.73 + 0.95 
1.44 + 0.74 

8.29 + 7.19 
2.03 + 1.29 
7.16 + 1.1 
0.37 + 0.49 
23.4 + 26.9 
64.0 + 15.3 
4.11 + 3.27 
4.66 + 5.68 
12.6 + 9.9 

6.49 + 2.24 
8.23 + 3.18 
0.87 + 0.24 
5.54 + 2.87 
7.01 + 2.19 
9.82 + 7.75 

10.5 + 7.8 
17.3 + 14.6 

27.6 + 19.4 
0.87 + 0.4 



Note. — Results of the Power Ratio analysis fSection |4.4|) 
of the clusters. 



TABLE 6 
Power ratio results 



Sample 


N.o. objects 


P2/Po(xlO-7) 


P3/Po(xlO-7) 


0.069 < z <0.1 


28 


71 + 3 


1.7 + 0.3 


O.K z< 0.3 


55 


89 + 5 


2.2 + 0.5 


z > 0.3 


18 


65 + 7 


7.0 + 2.8 


cooling clusters 


31 


39 + 2 


0.29 + 0.07 



Note. — Average values of the Power Ratios derived in different 
redshift bins. 



containing the clusters with the highest luminosity sub- 
structure (high P3/P0, see Section [5.2. ip shows signifi- 
cant evolution with z. Slightly higher than that of the 



non-cooling sample. 

It is possible that this has to do with the incomplete- 
ness of the sample. A cut on P3/P0 selects a larger frac- 
tion of high- 2; clusters, since P3/P0 has a weak positive 
trend with z. High redshift clusters tend to be more lu- 
minous by selection and this could cause a small bias in 
the observed evolution. We find that our results for the 
clusters without cool cores are in good agreement with 
theoretical models of cluster evolution that are based on 
balancing gas cooling with feedback at the e ntropy set by 
the cooling threshold (|Voit et al.ll2QQl : IVoitll2QQ5h . This 
threshold Kc{T^t) is the entropy at which constant- en- 
tropy gas at temperature T radiates an energy equivalent 
to its thermal energy in time t. These models are also 
successful in explaining the observed Lx oc behav- 
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Fig. 9. — Power ratios P2/PQ vs redshift for all clusters. Low luminosity clusters are shown as red circles, intermediate luminosity as 
green stars and high luminosity as blue squares, cf. Figure [Jl Simulated cool core clusters are shown as filled gray circles connected by 
black lines. The simulated two-subcluster model is shown as purple triangles connected by black lines. Letters from A to O denote a few 
selected clusters (see Figure [Jl). 




Fig. 10. — Power ratios P3/P0 vs redshift for all clusters. Low luminosity clusters are shown as red circles, intermediate luminosity as 
green stars and high luminosity as blue squares, cf. Figure [Jl Simulated cool core clusters are shown as filled gray circles connected by 
black lines. The simulated two-subcluster model is shown as yellow triangles connected by black lines. Letters from A to O denote a few 
selected clusters (see Figure [Jl). 



ior for clusters, contrary to early assumptions of pure 
gravitational self-similar collapse, where Lx cx was 
originally expected. 

5.3.1. Temperature correlation 

Since the temperature two-point correlation function is 
new, we illustrate its power for four representative clus- 
ters that clearly show specific characteristics. A{rk) is 
calculated for all spatial points, with 10 kpc resolution, 
within radius of 500 kpc of the cluster centroid. 

First, we consider a well-known cooling core cluster, 
Abell 1835 (see Fig. [HI first panel). Here, A{rk) is 



shown as a function of distance scale, r. As expected, 
there is considerable power at small spatial scales, im- 
plying a compact structure with temperature clearly dif- 
ferent from the rest of the cluster. There is a sharp 
drop after 500 kpc where the core is no longer visible. 
Another extreme example is MS1054.4-0321 (Fig. [131 
second panel), a highly substructured cluster consisting 
of multiple components. Here, the amplitude is high, 
distributed more uniformly over a large range of spatial 
scales, suggesting multiple high-luminosity components 
with different temperatures. As a third example, we 
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Fig. 11. — P3/P0 values for L^oi > 4 10^^ erg s ^ clusters up to 2; = 0.6 averaged in 5 redshift bins. Best fit line (solid line) is shown 
along with interpolated trends from the simulated clusters (dashed lines, red dashed line has lowest x^)- Simulated cool core clusters are 
shown as filled gray circles and simulated irregular (3-core) clusters are shown as yellow stars. 



TABLE 7 
Power ratio comparison 



Name 


P2/P0(X10-^) ^ P2/P0(X10-^) 


P3/Po(xlO-^) ^ 


P3/Po(xlO-7) 


A1413 


58.0 


192 ±3 


0.0525 


0.19 ±0.1 


A2034 


14.1 


84.5 ±8.4 


0.476 


0.41 ±0.21 


A2218 


27.4 


72.8 ±5.6 


0.233 


1.14 ±0.36 


A1914 


15.7 


32.4 ± 1.0 


0.709 


2.45 ± 0.26 


A665 


19.4 


44.8 ±3.5 


2.43 


7.75 ± 0.77 


A520 


40.7 


65.3 ± 11.2 


2.43 


3.87 ±0.46 


A963 


5.03 


8.14 ± 1.32 


0.342 


1.34 ±0.67 


A773 


46.9 


83.2 ± 11.2 


-0.125 


0.82 ± 0.39 


A2261 


4.57 


12.9 ± 11.4 


0.201 


2.79 ±1.7 


A2390 


58.0 


149 ±7 


0.291 


2.97 ±0.86 


A267 


61.4 


95.8 ±9.2 


-0.306 


0.92 ± 0.64 


RXJ2129.6+0005 


17.6 


66.1 ± 1.4 


-0.0814 


0.21 ±0.03 


A1758 


188.0 


502 ± 13 


1.06 


0.78 ± 0.32 


ZW3146 


4.42 


11.6 ± 1.3 


0.078 


0.37 ± 0.08 


MS2137-23 


1.81 


1.47 ±0.96 


0.00772 


0.37 ± 0.49 


CL0016+16 


46.4 


114 ± 17 


0.316 


5.54 ±2.87 


MS0451. 6-0305 


66.1 


82.3 ±5.1 


2.19 


7.01 ±2.19 


MS1137.5+6625 


5.24 


11.5 ± 11.9 


0.115 


10.5 ±7.8 


MS1054.4-0321 


150.0 


220 ± 73 


10.3 


17.3 ± 14.6 


WARPJ0152.7-1357 


264.0 


2880 ± 260 


12.4 


27.6 ± 19.4 


CLJ1226.9+3332 


-0.821 


11.1 ±3.2 


0.77 


0.87 ±0.4 



Note. — Co mparison of the Power Ratios in our analysis to the results using Chandra data 
bv lJeltema~a l. (2005) 

^ Derived bv lTeltema et al.l (|2005l ) from Chandra data. 



also show the correlation function for a nearly isother- 
mal cluster, A1689 (Fig. [131 third panel), where small 
temperature fluctuations lead to lower values across the 
range of spatial scales. Finally we show the correlation 
function for A2142 (Fig. [131 fourth panel), where an 
offset of the cluster core causes the sharp drop seen in 
A1835 to be absent. The drop is gradual, implying an 
overall asymmetry in the temperature structure of the 
cluster. 

The temperature correlation function A^r^) is inte- 
grated over 250-500 kpc {Ai) and over 750-1000 kpc {A2) 



and this is displayed in Figure [TH plotted against red- 
shift. This shows the magnitude of temperature gradi- 
ents over small and large scales respectively. These dis- 
tance ranges are chosen since for the small scales, 250- 
500 kpc, the cooling cores dominate the statistic. To 
distinguish these fluctuations from larger scale, merger 
related disturbances we compare it to the 750-1000 kpc 
range where the core cannot be included (since we use a 
500 kpc aperture). 

Strong cooling clusters such as A1835 and A2204 show 
large values of Ai since the cool cores dominate the small 
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Fig. 12. — Luminosity temperature relation for low redshift clusters (z < 0.2, red circles), high redshift clusters (z > 0.2, green stars) 
a nd cooling core c l uster s (blue squares). Best fit Lx-T relations are shown as solid lines and the non-cooling and cooling samples at z = 
of I Allen &: FabianI (|1998l ) are shown as dashed lines. 



TABLE 8 

Luminosity-temperature relation 



Sample 


L6 


aLT 


asyst(\og(T)) 


0.069 < z <0.2 


6.5 ±0.4 


2.98 ±0.19 


0.046 


z > 0.2 


8.9 ±0.7 


2.62 ±0.21 


0.050 


cooling clusters 


20.3 ±2.0 


2.69 ±0.18 


0.075 



Note. — Best fit parameters for the Lx-T relation. In addi- 
tion to the measurement errors, a systematic scatter in log(T) 
has been added in order to achieve /dof. = 1. The scatter in 
temperature in the non-cooling sample corresponds to ~ 0.7 
keV at 6 keV. 



scales. The same objects show low values of A2 since 
these clusters are largely isothermal when the core is ex- 
cluded. The same is true for our simulated cool core clus- 
ters, while the small scales cannot quite be resolved at 
high z. In contrast, disturbed clusters such as RXJ0658- 
55 or MS 1054.4-0321 show generally lower values of Ai 
and larger values of A2 whereas nearly isothermal clus- 
ters such as A1689 or A1413 show low values for both. 
With the exception of a few low luminosity clusters at low 
redshift, large scale temperature structure {A2) can be 
seen to increase slightly with redshift above z = 0.2 when 
compared to the simulated clusters. There is also an 
apparent drop in the temperature structure for smaller 
scales (Ai). 

We note that the decrease in Ai with redshift can be 
partly due to smoothing of the core caused by the PSF 
and loss of photons as seen in the evolution in the simu- 
lated clusters in Figure [TH However, it is unlikely that 
this effect is the cause of the increase in A2 since temper- 



ature features are likely to be washed out by this effect. 
This is seen in the trend of A2 for the simulated clusters. 

In Figure [15] we plot A2 against Ai to show the sepa- 
ration of "relaxed" cool core clusters (lower right) from 
disturbed clusters (upper left). 

5.4. Notes on individual dusters 

In this section we describe breifiy a few selected clus- 
ters in our sample, comparing our findings with previous 
results. The luminosity and temperature maps of these 
clusters are presented in Figures [16] to [32] 

5.4.1. Abell 1068 

A high resolution temperature map of Abell 1068, 
using Chandra data, was previously presented in 
IWise et aTl (|2004l ) who find an large degree of substruc- 
turing within the central 80 " x 80 " with temperatures 
varying between 2.5 and 6 keV. We find that this agrees 
well with our map in Figure [161 although the spatial res- 
olution is not as high as in the Chandra map. The size of 
the 80 " X 80 " region can be estimated using 4 times the 
PSF width in Figure [16] In a radial te mperature pro- 
file derived from the XMM- Newton data, ISnowden et al.l 
(2008) find a central temperature of 3.1 keV increasing to 
5.3 keV around 0.4 rsoo and dropping to 2.2 keV around 
0.8 rsoo (~ 800 kpc). In our map we find no indication 
of such a drop in T. We note, however, that due to the 
low cluster flux at these radii, the error on the temper- 
ature in the map is large (~ 3 keV). There may also be 
some associated systematic effects (see Section [6]). The 
average temperature of the cluster that we find from a 
simple spectral analysis is 3.81^o 07 keV. 
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TABLE 9 

Redshift dependent luminosity-temperature relation 



Sample Lq aLT (3lt crsyst(log(T)) 



non-cooling 6.8 ± 0.5 2.87 ±0.14 0.50 ± 0.34 0.037 

cooling clusters 18.8 ± 5.9 2.79 ± 0.26 0.62 ± 1.47 0.067 

high Lboi NC 7.6 ±0.7 2.81 ± 0.26 0.48 ± 0.32 0.029 

high P3/P0 NC 6.4 ± 0.6 3.02 ± 0.18 0.82 ± 0.42 0.035 



Note. — Best fit parameters for the Lx-T relation in Eq. 1101 In addition 
to the measurement errors, a systematic scatter in log(T) has been added in 
order to achieve /dof. — 1. The scatter in temperature in the non-cooling 
sample corresponds to ~ 0.5 keV at 6 keV. 
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Fig. 13. The temperature correlation function, A{rk) plotted for clusters A1835 (top left), MS1054.4-0321 (top right), A1689 (bottom 
left) and A2142 (bottom right). 



5.4.2. Abell 1413 

Early studies of the temperature profile of Abell 1413 
using XMM- Newton data find a central temperature of 
appr oximately 7 keV d eclinin g to 6 keV out t o 1 Mp c ra- 
dius (jPratt Arnaudll2QQ2f ). ISnowden et~an (|2QQ8h find 
a similar trend with a central temperature of 7.7 keV 
declining to 5.7 keV out to ^ 1 Mpc. Our temperature 
map (Figure [T7|) exhibits a large degree of small scale 
substructure in this otherwise smooth surface brightness 
distribution with temperatures varying between 7 and 
9 keV. We do not see a clear temperature drop with 



increasing radius but rather a large degree of asym- 
metry. We find an average spectral temperature of 
7.30t^;{^ keV. 

5.4.3. Abell 1689 

The temperature map of Abeh 1689 (Figure [18]) 
is in good agreement with our earlier analysis of 
the XMM-Newton da ta (lAndersson fc M adeiskl l2QQ4l : 
lAndersson et al.1 I2QQ7D . ISnowden et all ((20081) find a 
slightly higher temperature of 11.6 keV around 0.1 rsoo 
(~ 200 kpc) whereas at other radii the results are consis- 
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Fig. 14. — Temperature correlations for scales 250 — 500 kpc (Ai) vs redshift (left panel) and 750 — 1000 kpc (A2) vs z (right panel). 
The different symbols represent luminosity (red circles > 10^^ erg s"-*^, green stars > 4 10^^ erg s"-*^, and blue squares > 2 10^^ erg s-i). 
Simulated cool core clusters are shown as filled gray circles. Letters from A to O denote a few selected clusters (see Figure [7}. 




Fig. 15. — Temperature correlations for scales 750 — 1000 kpc (A2) plotted against 250 — 500 kpc (Ai) scales. Here the different symbols 
represent redshift (red circles 0.069 < 2; < 0.1, green stars 0.1 < 2; < 0.3, and blue squares z > 0.3). Simulated cool core clusters are shown 
as filled gray circles. Letters from A to O denote a few selected clusters (see Figure [7}. 

tent with ours with fluctuations from 8 to 10 keV. The a simple spectral fit we find the average temperature of 
average spectral temperature is 9.07to;l2 keV. A1750 to be 4.451q;i3 k^"^- 



5.4.4. Abell 1750 

The XMM-Newton data for Abe ll 1750, consis t ing o f 
two sub clusters, was analyzed by iBelsole et alJ (|2004l ) 
who find temperatures varying between 3 and 6 keV 
within the double-cluster. The temperature structure 
overall is in good agreement with our temperature map 
while we find temperatures varying from 4.5 to 8.5 keV. 
This is likely related to the bias described in Section [6l4l 
We find that the larger subcluster (pictured center in Fig. 
[T9|) is slightly hotter (5.5 keV) compared to the north- 
ern subcluster (4.5 keV) with the gas in betw een having 
a temperature of 8 keV. iBelsole et al.l (|2004| ) find these 
temperatures to be 2.8, 3.9 and 5.1 keV respectively. In 



5.4.5. Abell 1758 

Abell 1758 is a complicated system consisting of a 
northern and a southern component, separated by 2 Mpc 
in the plane of the sky, each of which, in turn, consists of 
two separate subclusters. We have limited our analysis 
to the northern component A1758N (See Fig. [2Q|) . The 
Chandra and XM M-Newton data of AI 758 have been 
analyzed earlier by iDavid fc Kempnerl (|2004l ) who con- 
clude that A1758N is in the late stages of a large im- 
pact parameter merger of t wo 7 keV clusters. The ha rd- 
ness ratio map presented in iDavid fc Kempnerl (|2004l ) is 
qualitatively in good agreement with our temperature 
map. For their spectral fits in four separate regions they 
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find temperatures of 7.2, 6.6, 7.2 and 9.8 keV for the 
"northwest wake" , "southeast wake" , "core" and "halo" 
using the XMM-Newton data. The values for the same 
regions in our temperature map are systematically ap- 
proximately 1 keV higher. In our overall spectral fit we 
find a temperature of 9.16lo.43 

5.4.6. Abell 1775 

Abell 1775 is a nearby {z = 0.072), low mass clus- 
ter with high velocity dispersion, containing a pair of 
massive galaxies at the center (Figure [21]). Studying 
the velocity distribu tion of the galaxies in this cluster, 
lOegerle et al.l (|l995f ) conclude that it, in fact, consists of 
two smaller subclusters that are in the process of merg- 
ing. In our temperature map, we can see an arc of the 
coldest gas extending to the east of the cluster center, 
possibly a remnant of a subcluster core. In this low mass 
cluster, the surface brightness becomes low outside 500 
kpc and the temperature map outside of this radius is 
not reliable. Within 500 kpc we find the gas temper- 
ature varying in the 4 to 7 keV range, biased upward 
by effects described in Section 16. 4[ Adapting the prior 
range of the temperature distribution in this case to bet- 
ter accommodate for low temperature easwould give a 
more accurate result. ISnowden et al.l (|2008l ) find a cen- 
tral temperature of 4.3 keV decreasing to 3.5 keV around 
0.4 rsoo 600 kpc). In our overall spectral fit we find 
a temperature of 3.70lo!o6 keV. 

5.4.7. Abell 1835 

Abell 1835 is a massive cluster w ith a pronounced 
cool c ore and a relaxed morphology. iMajerowicz et all 
(|2002l ) find a central temperature of 4.1 keV increas- 
ing to 8.5 keV near 1 Mpc using XMM-Newton data. 
ISnowden et al.l (|2008l ) find a central temperature of 
6.1 keV, increasing to 10.5 keV just outside the clus- 
ter core and decreasing to 7.7 keV out to 1 Mpc using 
the same data. This is in good agreement with our tem- 
perature map shown in Figure [22l In a simple spec- 
tral fit we find the average temperature of A1835 to be 



5.4.8. Ahell 2029 

Using an early Chandra observation of A2029, 
iLewis etliD (|2002l ) find a central temperature as low as 
4 keV within 10 kpc increasing to a pproximately 9 keV 
out to 300 kpc. Using the same data,IClarkeelalJ (|2004f ) 
find a disturbed core within a 20 kpc radius with some 
spatial correlation to the central radio emission. In our 
temperature map in Figure [23] we find no significant cen- 
tral substructuring o n the scales we are probing w ith the 
XMM-Newton data. iBourdin fc Mazzottal (|2008[ ) find a 
central temperature of 5.7 keV and a large degree of tem- 
perature substructure with fluctuations between 7 and 
9.5 keV in their temperature map out to 600 kpc, gener- 
ated using the XMM-Newton data. While in good quan- 
titative agreement with the radial trend in our tempera- 
ture map, we find no sign of the hot or cold substructures 
in our map. In our overall spectral fit we find a temper- 
ature of 6.88t^;J^ keV. 

5.4.9. Ahell 2142 



iMarkevitch et all (|2000[ ) first argued that A2142 con- 
tains a "cold front" , a steep nonaxisymmetric X-ray flux 
gradient and a steep temperature gradient originally be- 
lieved to be due to a low-entropy remnant of a r ecently 
merged subcluster. iTittlev fc HenriksenI (|2005l ) argue 
that this "cold front" is due to oscillatory motions of 
the core after a merger event. Our temperature map 
(Figure [24]) shows a central temperature around 7.5 keV 
with a plume-like structure of 8.5 keV gas extending 
north and a steep temperature gradient, aligned with 
the sharp trend in surface brightness to the north-west, 
increasing to 11.5 keV. T his agree well with the results of 
[Markevitch et al.[ (|2000l ). with the exception of the cen- 
tral 5 keV gas which we do not detect. We find an overall 
spectral temperature of A2142 to be 8.15^Q"2g keV. 

5.4.10. Ahell 2204 

Abell 2204 is a massive cool core cluster with a regular 
morphology, ganders et al. (2005) find the temperature 
of the c ore region to b e 3.26d =0.20 keV. The temperature 
map of [Sanders et al.[ (|2005[ ) shows a gradually increas- 
ing temperature to ~ 10 keV at 100 kpc, continuing up 
to 16 keV around 400 kpc after which it drops again to 
around 7 keV. The trend is continuous, with the excep- 
tion of an arc-shaped region of hot gas extending around 
the northern part of the cluster at 50 kpc where the tem- 
perature suddenly jumps to 14 keV. Our map in Figure 
[25l shows a similar trend with a central temperature of 
4.5 keV gradually increasing to 10 keV around 300 kpc 
where after it drops again to 7 keV at 1 Mpc. We find 
no sig ns of the hotter gas around 50 kpc. [Reiprich et al.[ 
(|2008[ ) also find a similar radial profile using data from 
Chandra, XMM-Newton and Suzaku telescopes with a 
sharp increase in temperature outside of the core and a 
gradual decline out to the virial radius where a temper- 
ature of 4 keV is measured. In our overall spectral fit we 
find a temperature of 6.44^Q Qg keV, dominated by the 
emission from the core. 

5.4.11. Ahell 520 

A detailed temperature map of A520, using a 
short Chandra observ ation, was first presented in 
[Markevitch et al.[ ([2003') . The map shows the southwest- 
ern cluster core having a temperature of approximately 
5 keV while the main body of the cluster exhibits an 
irregular morphology with temperatur es varying from 6 
to 14 keV. Using a longer observation, [Markevitch et all 
(|2005[ ) detect a temperature drop in the southwest quad- 
rant, from 11.513'^ ~ ^00 kpc distance from 
the center to 4.8^0 j keV at 600 - 1000 kpc. We do not 
find any evidence for such a temperature drop. This is 
most likely due to the fact that the surface brightness is 
low in this quadrant at this large radius. Directly south 
of the main cluster, we see a region of possibly shocked 
gas at 11 keV and this appears to be the source of the 
hottest emission in this cluster (see Figure [26]) . In the 
overall spectral fit we find this cluster has an average 
temperature of 8.45to:26 "^^^ • 

5.4.12. MSO451. 6-0305 

The temperature map in Figure [27| reveals a low- 
temperature (7 keV) structure east of the main cluster 
(10 keV). It is possible that this component represents a 
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subcluster in the process of merging with the main clus- 
ter, or perhaps a less massive cluster in projection. This 
feature can be seen in the maps of surface brightness 
presented in lDonahue et al] ([2003), where the soft-band 
image is seen to extend furt her east than the hard-band 
map. iDonahue et al.l (|2003l ) further find a central tem- 
perature (r < 31^') of 10.31^ 5 outer tem- 
perature (31'' < r < 85'') of 9.8ti j keV. In our overah 

spectral fit we find a temperature of 10.02to.6o 
good agreement with their results. 

5.4.13. MS10544-0321 

MS1054. 4-0321 has been s tudied prev i ously with 
Chandra and XMM-Newton. iGioia et all (|2004D find 
an overall temperature of 7.2^0.7 keV using the XMM 
data while iJee et al] (|2005f ) find an overall temperature 
of 8.9to^ keV using the Chandra data. iJee et all (|2005f ) 
also find the temperature of the western peak to be 
7.5t]^"2 keV, the central peak to be 10. TI^;^ keV and an 
overall positive gradient in their temperature map from 
west to east. We confirm the presence of a temperature 
gradient, seen rising from 9 to 11 keV in Figure[28l In our 
overall spectral fit we find a temperature of 9.20t{ o3 
in good agreement with the Chandra results. 

5.4.14. MS2137-23 

MS2 137-23 is a strong gravitational lensing cluster 
with a well defined central cool core. In the tempera- 
ture map in Figure [29] we see the temperature increase 
from 5 to 8 keV within 500 kpc radius. The map exhibits 
a low degree of asymmetry which appears to increase 
with radius. In a depro.jec ted temperature profile us- 
ing a Chandra observation, lArabadjis et al.l (|2Q04) find 
a central (~< 50 kpc) temperature of 4 keV rising to 
around 7 keV at 150 kpc radius and decreasing again at 
larger radii. An offset datapoint at 250 kpc radius shows 
a temperature of 12.5 ± 4.5 keV. Our overall spectral fit 
shows a temperature of 4.67^0. 19 keV, clearly showing 
the dominance of the cool core emission. 

5.4.15. RXJ0658-55 

For a detailed dis cussion on the temperat ure structure 
of RXJ0658-55, see lAndersson et al.l (|2007[ ). Our maps 
of luminosity and temperature are shown in Figure [30l 

5.4.16. RXJI347-II45 

The temperature map in Figure [31] shows a central 
temperature of 9 keV increasing to 11.5 keV within 500 
kpc and decreasing again towards 1000 kpc radius. The 
temperature map further shows a lack of the hottest 
emission in the north quadrant of the cluster. We com- 
pare o ur results with the radial temperature profile and 
map of lGitti fc Schindlerl (|2004f ). derived from the XMM- 
Newton data. The temperature profile shows an overall 
good agree ment with our results, ho wever, the tempera- 
ture map of lGitti fc Schindlerl (|2004[ ) identifies a region of 
hot (~ 17 keV) emission approximately 150 kpc south- 
east if the core. This region appears in our map as a 
region of 11.5 keV emission aligned with the elongation 
of the surface brightness to the SE. Our map also shows 
less substructure overall. We find this massive cluster 
has an overall spectral temperature of 11.44^0^29 keV. 



5.4.17. ZW3146 

The temperature map of the cool core cluster ZW3146 
(Figure [32]) reveals a central temperature of 5.5 keV in- 
creasing to approximately 9 keV at 500 kpc. At larger 
radii the temperature appears to increase further, how- 
ever, here the measurement i s unc ertain. Using XMM- 
Newton data, Kau sch et al.l (|2007l ) find a temperature 
gradient from 4.7 to 7 keV up to 1 Mpc with most of the 
colder emission being localized within 100 kpc. Their 
temperature map shows fluctuations between 4.5 and 9 
keV in within 200 kpc. We do not see these fluctuations 
and conclude that the cluster has a relatively relaxed 
morphology with a smooth temperature gradient. The 
surface brightness contours have tighter spacing to the 
south-east which may indicate recent core motion rela- 
tive to the cluster as a whole. Our overall spectral fit 
shows a temperature of 6.21 to. 10 keV. 

6. SYSTEMATIC EFFECTS 

6.1. Selection bias 

Since the cluster sample was selected using all the pub- 
lic XMM data available at the time, it does not consist of 
a representative subset of clusters in general in the uni- 
verse. However, many studies, and observation propos- 
als, are based on either the most regular known clusters 
or the most complex. Hence, the sample should contain 
both these extremes of the cluster population which is 
what is compared in this work. 

6.2. Background modeling 

The X-ray background is modeled as opposed to sub- 
tracted which is the usual procedure in standard spectral 
modeling of extended X-ray sources. This introduces an 
uncertainty different from the unsubtracted background 
flux which is common when dealing with combinations of 
several sources of background with different spatial and 
spectral distributions. Here, the distortions of different 
backgrounds are modeled correctly, but instead, there is 
the possibility of cluster particles mimicking the back- 
ground flux. This effect is small since the backgrounds 
are present across the entire field which is not the case 
for the majority of the cluster emission since all clusters 
are modeled out to at least 1 Mpc where the cluster emis- 
sion tends to be faint. The low cluster flux at this radius 
does, however, impact on the precision of the temper- 
ature measurements at large radii. The nature of this 
uncertainty is seen in some temperature maps of median 
temperature as it reverts to the mean of the prior, which 
in this case is 7.75 keV. 

The spectral shape of the instrumental background 
also means that it is not easily modeled by cluster parti- 
cles. The particle background has a nearly flat spectral 
shape across the entire energy range whereas the X-ray 
photons are modulated by the energy dependent effective 
area of the X-ray mirrors. This makes it statistically 
favorable to model these events using the background 
model which is not affected by mirrors as opposed to 
using the SPI particles. 

In order to estimate the effect of the background mod- 
eling on the temperature at large radii we compare the 
temperature correlation results, Ai and A2, for the sim- 
ulated regular clusters with and without background. 
These are shown as filled circles in dark and light gray 
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respectively in Figures [T4l and [T5l The inclusion of back- 
ground leads to an overestimation of approximately 25% 
in Ai and 35% in A2. This is comparable to the l-cr 
statistical errors and since the effect appears to be rel- 
atively independent of redshift we conclude that it does 
not affect the evolution of the temperature correlation. 

6.3. Point-sources 

Point sources can significantly alter maps of cluster lu- 
minosity and temperature. In the cases where we can 
identify the point sources from the original event file we 
remove the flux from these sources using the technique 
described above (Section [4.H Figure [5]). When a point 
source cannot be resolved it can have a large impact on 
the temperature maps even if the effect on the luminos- 
ity map is not visible. Ideally, Chandra data of the same 
cluster should be used in combination to identify the po- 
sitions of point sources and to verify the accuracy of tem- 
perature features. Chandra data is however not available 
for all clusters in the sample and this type of analysis is 
beyond the scope of this work. 

6.4. Temperature bias 

As mentioned earlier, the average temperature from 
the SPI modeling becomes slightly higher than when 
standard analysis routines are applied. This effect stems 
from the fact that we are using multiple particles in our 
modeling. Since the particles are extended and hundreds 
of particles are required to model a single cluster there 
will be multiple overlapping cluster particles at any given 
position. Each particle has a separate temperature and 
set of elemental abundances and this means that the 
model is inherently multi-phase. There are many com- 
binations of thermal spectra that can be constructed to 
agree with the observed data and the final model consists 
of a wide range of phases at each spatial position. This 
range is very much dependent on the prior range for the 
parameter in question. Here, these ranges are visualized 
by making use of the median of the range at each spa- 
tial coordinate and this value again depends on the prior 
range. Since the measured differences of thermal spec- 
tra in XMM become smaller with increasing temperature 
the median is naturally biased towards the higher part of 
the prior. This may lead to a potential luminosity bias in 
the bolometric correction since the luminosities are cal- 
culated separately for each smoothed particle which will 
cause a higher luminosity if the temperature is higher. 
We estimate this offset to be less than ^ 20 % based 
on a comparison when using the temperature from the 
"standard" analysis in the bolometric correction. 

6.5. Modeling artifacts 

Since the modeling itself is based on particles, even if 
between 100 and 1000 particles are required, there is a 
possibility for particle artifacts appearing in the lumi- 
nosity and temperature maps. The cluster particles are 
constantly moving within the parameter space and inter- 
changing positions with each other, and averaging over 
many Markov chain samples reduces this effect. How- 
ever, since we are limited by CPU resources (see also 
Section 14. 2p and only 250 samples are used in the er- 
ror calculation these artifacts may still show up in the 
maps as lumpy structures and this may have a small ef- 



fect increasing the power ratios. In comparing power ra- 
tios calculated from maps using 1500 samples with ones 
generated using 250 samples we establish that this effect 
corresponds to about a 20 % increase for ratios (P2 / or 
P3/P0) between 10 ~^ and 10~^ and less for larger values. 
However, this is not sufficient to explain t he discrepancy 
in the comparison with the values from iJeltema et al.l 

dlooll). 

6.6. Trends of the simulated clusters 

In many of the plots, the trends in the relevant param- 
eters with z for the simulated cool core clusters appear 
similar to the evolution for the sample as a whole. In 
most cases, the trend for the simulations only represents 
the minimum disturbance that can be measured and the 
similarity does not necessarily mean that the evolution 
is not real. In fact, we see from our simulations of ir- 
regular clusters that the trend for a cluster with high 
P2/P0 or P3/P0 is much more stable than that of the 
circularly symmetric clusters. Ideally, to make a more 
rigorous claim of evolution, clusters of many different 
morphologies should be simulated and the apparent evo- 
lution analyzed separately. This will be the subject of 
future work and we are limited here, again, by CPU con- 
straints. Markov chains should also be run for several 
more iterations to limit model noise. 

7. DISCUSSION AND CONCLUSIONS 

We have applied a novel technique for X-ray galaxy 
cluster modeling to a large sample of XMM observed 
galaxy clusters in order to identify clusters with cooling 
cores and to study the redshift dependence of cluster sub- 
structure and the luminosity-temperature relation. We 
also developed a new statistic to quantify thermal sub- 
structure and to identify and distinguish cooling core 
clusters from isothermal and thermally disturbed clus- 
ters. We note here that while this technique is impor- 
tant in identifying clusters for which mass measurement 
can be used for cosmological applications, in its present 
form it is not capable of estimating cluster masses. How- 
ever, a straightforward extension of it, with additional 
assumptions, should allow the determination of three- 
dimensional density and temperature of the cluster, and 
thus the mass of both gas and dark matter content. 

The number of cooling core clusters appears to be lower 
at high z with 32 % identified as cooling clusters at z < 
0.1, 36 % at 0.1 < z < 0.3 and 11 % at z > 0.3. This 
suggests that the formation of cooling cores in clusters is 
a recent phenomenon. However, due to the low number 
of clusters and various selection effects these percentages 
are uncertain. 

7.1. L-T relation 

The derived luminosity-temperature relation for clus- 
ters clearly shows a large offset for cooling core clusters, 
as expected. This highlights the importance of apply- 
ing a redshift dependent criterion when excluding these 
sources. In contrast to previous work we do not attempt 
to correct for cooling cores by applying some standard- 
ized central profile. Inst ead we fit the relation for cooling 
clusters separately (cf. lAllen fc FabianI Il998l ) and find 
a similar slope of ~ 3 but with a normalization ~ 3 
times higher than for the non-cooling sample. The nor- 
malization of our Lx-T relation for cool core clusters 
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is ~ 5 0% higher for our sam ple compared to the find- 
ings of I Allen FabianI (|1998[ ). This is likely due to the 
use of ROSAT and Asca data to determine temperatures 
and luminosities, where cool core clusters would hard to 
model with the low spatial resolution and the limited 
energy band. 

We also note that the scatter in the Lx-T relation is 
greater for cooling core clusters. We detect weak evolu- 
tion in the non-cooling Lx — T relation oc (\j^z)^-^^^^-^^ 
in good agreement with self similar predictions and with 
an average slope of a^T = 2 .87 ± Q-l^' ^^} ^ con- 
trast to the findings of Vikhli nin et al.l (|2QQ2l ) who detect 
strong evolution oc (1 + Due to the limited ob- 

servable volume available at low high luminosity clus- 
ters are observed at mostly high redshifts. This could 
create a aLT-^LT bias. We note also that when study- 
ing Plt for the cluster sample with Ps/Pq > 10~^, the 
sample shows significant evolution. This may be an im- 
portant effect to account for in large surveys relying on 
the Lx-T relation to measure the mass function of clus- 
ters over a range in z. Clearly, the Lx-T relation has 
to be measured in greater detail with larger samples of 
well exposed clusters in order to reach a consensus of its 
evolution. To measure the mass function of clusters to 
constrain cosmological parameters the knowledge of this 
evolution is crucial. 

7.2. Structure quantification 

The power ratio analysis was applied uniformly for the 
101 clusters within a 500 kpc radius. We find no sig- 
nificant increase in P2/P0 with z but detect an increase 
in P3/P0 implying a higher level of substructure at high 
redshift. Our sample of cooling core clusters shows sig- 
nificantly lower values of P2 / Pq and P3 / Pq than the non- 
cooling sample, which is not surprising since cooling clus- 
ters generally are the most dynamically relaxed clusters. 

Comparing this quantification of cluster structure with 
simulations of spherical clusters we show that all clusters 
in our sample are far from spherical, as determined by 
the value of P2/P0 and should not be modeled as such. 
This is easily realized when viewing the luminosity maps 
in Figures [T6II32I which exhibit a large degree of ellipticity. 
The large amount of features visible in the temperature 
maps confirm that a large number of processes govern 
the properties of the intra-cluster medium, that cluster 
merging can give rise to highly complex thermal struc- 
ture and that temperature seldom is a simple function of 
radius. 

The SPI method inherently lends itself to a study of 
temperature structure of clusters, and to take advantage 
of it, we design a statistic, the "temperature two-point 
correlation", capable of distinguishing between cooling 
clusters, isothermal clusters and clusters with signifi- 
cant temperature differences over large distances, such as 
merging clusters. The analysis using this function hints 
at an increase in the large scale thermal structure of clus- 
ters with redshift, indicating more frequent merging in a 
denser universe. This indication could be confirmed us- 
ing a joint approach utilizing both XMM- Newton and 
Chandra datasets. 
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Fig. 16. — 2x2 Mpc images of the bolometric luminosity, in units of 10^^ erg s"-*^ {")~'^ (left panel), and temperature, in keV (right 
panel), for Abell 1068. The size of the 500 kpc region used for power ratio analysis and two-point temperature correlation analysis is shown 
along with the size of the PSF (10^^ radius) . 
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Fig. 17. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 1413. 
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Fig. 18. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 1689. 



Fig. 19. 
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2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 1750. 
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Fig. 20. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 1758. 
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Fig. 21. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 1775. 
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Fig. 22. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 1835. 
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Fig. 23. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 2029. 
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Fig. 24. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 2142. 
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2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 2204. 
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Fig. 26. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of Abell 520. 
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Fig. 27. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of MS0451. 6-0305. 
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Fig. 28. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of MS1054. 4-0321. 
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Fig. 29. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of MS2137-23. 



5 


1 






1 






0.0005 



0.002 0.0035 3.5 9 9.5 10 10.5 



Fig. 30. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of RXJ0658-55. 
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Fig. 32. — 2 Mpc x2 Mpc field showing luminosity (left) and temperature (right) maps of ZW3146. 
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